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Using large-scale dynamical cluster quantum Monte Carlo simulation, we study the Lifshitz tran- 
sition of the two dimensional Hubbard model with next-nearest-neighbor hopping (£'), chemical 
potential and temperature as control parameters. At t' < 0, we identify a line of Lifshitz transition 
points associated with a change of the Fermi surface topology at zero temperature. In the overdoped 
region, the Fermi surface is complete and electron- like; across the Lifshitz transition, the Fermi sur- 
face becomes hole-like and develops a pseudogap. At (or very close to) the Lifshitz transition points, 
a van Hove singularity in the density of states crosses the Fermi level. The van Hove singularity 
occurs at finite doping due to correlation effects, and becomes more singular when t' becomes more 
negative. The resulting temperature dependence on quantities such as the bare d-wave pairing sus- 
ceptibility close to the Lifshitz points is significantly different from that found in the traditional 
van Hove scenarios. Such unambiguous numerical observation of the Lifshitz transition at t' < 
extends our understanding of quantum critical region of the phase diagram, and shines lights on 
future investigations of the nature of the quantum critical point in the two dimensional Hubbard 
model. 

PACS numbers: 74.40.Kb, 71.10.Fd, 74.72.-h, 71.10.Hf 



I. INTRODUCTION 

The physical properties of high-T c cuprate supercon- 
ductors are extremely sensitive to doping. Experiments, 
such as the angle-resolved photoemission spectroscopy 
(ARPES) and quantum oscillation measurements 2 -, have 
clearly demonstrated a Fermi surface reconstruction as 
the doping concentration is varied. In overdoped mate- 
rials, ARPES and quantum oscillation studies have con- 
firmed the existence of a large Fermi surface compatible 
with band theory. As the doping is reduced, photoemis- 
sion data indicate that the Fermi surface breaks up into 
"Fermi arcs'—, and more recent studies further indicate 
that the "Fermi arcs" are actually part of closed hole 
pockets^ 5 -. The formation of small Fermi pockets, as the 
doping is reduced from overdoped to underdoped region, 
has also been substantiated by quantum oscillation mea- 
surements^. 

From a theoretical point of view, the change of the 
Fermi surface topology from a large surface to small 
"arcs" or "pockets", resembles a Lifshitz transition. Al- 
though proposed by Lifshitz in noninteracting Fermion 
systems decades age 8 , the Lifshitz transition has only 
recently begun to be considered as a quantum phase 
transition in strongly correlated electron systems^ — . It 
has been used to explain experimental data in high-T c 
cuprate superconductor s 15 ! 16 or heavy fermion systems^. 

Using large-scale dynamical cluster quantum Monte 
Carlo simulations 18 , a series of recent numerical 
works^r— mapped out the phase diagram of the hole- 
doped two dimensional Hubbard model near the quan- 
tum critical filling. Particularly, the superconducting 
dome in the proximity to the quantum critical doping 
has been identified 22 . At positive t' , there is a first-order 



phase separation transition occurring at finite temper- 
ature. The two phases being separated are an incom- 
pressible Mott liquid and a compressible Mott gas; these 
two phases are adiabatically connected to the pseudogap 
and the Fermi liquid states at t' = 0. The first-order 
line of coexistence terminates at a second order point 
where the charge susceptibility diverge s 21 ! 23 . As t' —> 0, 
this critical point goes continuously to zero temperature 
and thus becomes the quantum critical point (QCP) un- 
derneath the superconducting dome. Above the QCP, 
a V-shaped quantum critical region, characterized as a 
marginal Fermi liqui d 25 ! 26 with linear resistivity, sepa- 
rates the pseudogap and the Fermi liquid phases. Fur- 
thermore, when the next-nearest-neighbor hopping be- 
comes negative, t' < 0, there is indication of the Fermi 
surface topology change at zero temperature, and the fill- 
ing at which such changes occurs is an extension of the 
QCP at t' = to negative t'^i. 

Other numerical works also found similar feature of the 
Fermi surface reconstruction at negative t' by varying 
the hole doping concentration from overdoped towards 
half-filling. These include results from both dynamical 
cluster approximation by Gull et al*2£ — and results from 
the cellular dynamical mean-field theor)*2I by Sakai et 
al^-^i. 

Following such evidence, and in light of viewing 
the Fermi surface reconstruction as a Lifshitz transi- 
tion^ n i 13 i 15 i 16 , we performed systematical numeric stud- 
ies on the quantum critical phase diagram of the two 
dimensional Hubbard model at various t' with special 
attention to the region of t' < 0, which is relevant for 
the hole-doped cuprates. We find at t' /t < 0, as the 
doping concentration varies from overdoped to under- 
doped regimes, the Fermi surface changes its topology 
from electron-like with complete Fermi surface to hole- 
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like with pseudogap at the anti-nodal direction. Such a 
topological transition in the Fermi surface is a Lifshitz 
transition. It is furthermore associated with a van Hove 
singularity in the density of states crossing the Fermi 
level, at a doping which occurs very close to (if not at) 
the quantum critical point. 

The van Hove singularity crossing and the concomitant 
quantum critical point occur at finite dopings due to cor- 
relation effects even when t' = 0. Interestingly, we find 
the quantum critical phenomena prevail for negative t 1 
and the van Hove singularity defines a line of quantum 
critical points which extends from the QCP at t' = to 
t' < at higher dopings. The temperature dependence 
of these correlation effects close to the van Hove singu- 
larities, and its influence on quantities like the pairing 
polarization are very different from those found in the 
traditional van Hove scenario s 35 ! 36 . The QCP, and the 
van Hove singularity, have great impact on the conditions 
for pairing. At t' < 0, we find an enhanced tempera- 
ture dependence of the bare d-wave pairing susceptibility 
above the QCP that cannot be captured by the conven- 
tional BCS logarithmic divergence or the log-squared di- 
vergence found at a static van Hove singularity 35 . 

This paper is organized as follows. Section II outlines 
the model and the methods used in this study: the dy- 
namic cluster approximation (DC A) with weak-coupling 
continuous-time quantum Monte Carlo (CTQMC) as 
cluster solver. Section III and IV contain our numeri- 
cal results and discussion, beginning with the spectral 
function and dispersion at various doping concentrations 
and t'/t, and followed by a detailed account of other 
single-particle properties such as density of states and 
quasiparticle fraction, across the Lifshitz transition. We 
then provide results on the unique temperature depen- 
dence of the bare d-wave pairing susceptibility. We use a 
schematic quantum critical phase diagram of the model 
which summarizes our numerical results. We end with 
conclusions and an overview of open questions in Section 
V. 



II. FORMALISM 

In this work, we look for direct evidence of a Lifshitz 
transition in the spectral function of the two-dimensional 
Hubbard model 



(i) 



where c ]^ a ( c ka) 1S ^ e crea tio n (annihilation) operator 



for electrons with wavevector k and spin a, n. 
is the number operator, and the bare dispersion is given 

by 

e^, = — 2t (cos k x + cos k y ) — At' (cos k x cos k y — 1) , (2) 

with t and t' being the hopping amplitude between near- 
est and the next-nearest-neighbor sites respectively, and 
U the on-site Coulomb repulsion. 



We employ the DC A 37 i 38 with weak-coupling 
CTQMC 3 - 9 as a cluster solver. The DCA is a cluster 
mean-field theory that maps the original lattice onto 
a periodic cluster of size N c = (D is the dimen- 
sionality) embedded in a self-consistently determined 
host. It treats the spatial short-ranged correlations (up 
to L c inside a cluster) explicitly while approximating 
the long-ranged correlations with a mean-field. In this 
work we choose a square cluster with N c = 16, and 
the six independent momentum patches are centered 
at T = (0,0), M = (7r,7r), X = (tt,0), (7r/2,7r/2), 
(7r,7r/2), and (tt/2,0). We set the energy scale to 
At = 1, choose the interaction strength as U = 6t and 
study inverse temperatures up to f3 = 58/4t. The 
temporal correlations, essential for quantum criticality, 
are treated explicitly by the weak-coupling CTQMC 
solver for all cluster sizes. The solver expands in the 
Coulomb interaction diagrammatically and samples 
in time continuously. Different from the Hirsch-Fye 
algorith m 40 ! 41 it therefore has no Trotter error. 

We obtain high-quality estimates of the cluster 
self-energy Y,(K,oj) by employing the maximum en- 
tropy analytical continuation 42 (MEM) directly to the 
Matsubara-frequency self energies calculated from the 
DCA-CTQMC 2 !^^. The idea of performing the MEM 
on the self-energy is to normalize the non-Hartree part of 
the H(K,icj n ) by U 2 Xa,a, where x*,* = (n a n a ) - (n a ) 2 = 
n a (l — n a ) and is the local polarization of a single spin 
species a. The normalized spectrum of self-energy acts 
as a probability distribution: 



E(K,^ n )-E^ = f d ^ a(K,uj) 



(3) 



1. 



where a(K,uj) = E (K,u)/U Xa,a^ J duoa(K \uS) = 

7T 

1, using Xa,a obtained from the Monte Carlo process. 

To obtain the lattice self energy, E(/c,cj), we interpo- 
late the cluster self energy, Y,(K,oj). From the lattice 
self energy we can get the lattice single-particle spectral 
function, A(k, uo). An alternative way has been suggested 
by Stanescu et al^ 2 - 34 i 45 , where the cluster cumulant, 
M(K,u) = 1/(cj + fi — E(X,cj)), has been used for the 
interpolation to lattice quantities. Although the spec- 
tral functions produced by the two interpolation schemes 
have differences in the underdoped region, we found the 
two methods give the same quantum critical filling nj&. 



III. RESULTS 

Figure [1] shows the k-resolved spectral function eval- 
uated at Fermi level, A(k,uo = 0), for t'/t = and 
t'/t = —0.1. At t' = we observe a change of the 
Fermi surface topology from electron-like to hole-like as 
the filling increases from n = 0.75 towards half- filling. 
Such a change of topology of the Fermi surface is a Lif- 
shitz transition. It occurs at the same filling where the 
van Hove singularity crosses the Fermi leve l 36 ! 47 . The 
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FIG. 1. (color online) Zero frequency spectra A(k,ou = 0) for t' /t — (left) and t'/t — —0.1 (right) as a function of filling 
n. (Left panel) As n increases towards half- filling, the system undergoes a Lifshitz transition where the Fermi surface changes 
topology from electron-like to hole-like. The Lifshitz transition occurs at n c ~ 0.88, where the van Hove singularity crosses the 
Fermi level (as shown in Figure [2}. (Right panel) the Lifshitz transition occurs at n c ~ 0.835. Note that when t'/t = —0.1, the 
Lifshitz transition happens earlier than the t' = case. In either case, asn^ 1, the system enters the pseudogap phase with 
vanishing quasiparticle weight. 



corresponding filling is the critical filling, n c , of the Lif- 
shitz transition. In the case t' = we obtain n c « 0.88. 
At t'/t = —0.1 the Fermi surface changes from electron- 
like to hole-like and then disappears as n varies from 
0.8 towards half filling. The Lifshitz transition occurs 
at n c « 0.835, following the same criterion. Finally, for 
t'/t = -0.2, we find n c w 0.77 (not shown). For all 
t'/t < we studied, the system enters the pseudogap 
phase once the filling is larger than the corresponding 
n c . In the case of n = 0.92 for the left panel of Fig. [TJ 
and n = 0.87 for the right panel of Fig. [TJ the collapses 
of the single-particle spectral weight along the antinodal 
direction can be clearly seen. 

Figure [2] shows density of states iort'/t = (left panel) 
and t'/t = —0.1 (right panel), with the same set of fill- 
ings in Fig.[TJ For t' = 0, the van Hove singularity shifts 
from positive to negative frequency as the filling increases 
from 0.75 towards 1.0. The van Hove singularity is lo- 
cated at the Fermi level for the quantum critical filling 
n c w 0.88. These results are consistent with our previ- 
ous observations^*^. After passing through the critical 
filling, the weight between the lower Hubbard band and 
the upper Hubbard band is strongly suppressed, i.e. a 
pseudogap is formed. This parameter regime can thus be 
identified as the pseudogap region in the phase diagram. 
For t'/t = —0.1, the results in the right panel of Fig. [2] 
show a quantum critical filling n c ~ 0.835. Note, the van 
Hove singularity in single-particle density of states has 
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FIG. 2. (color online) The single-particle density of states for 
the parameters shown in Figure [1] (Left panel) t'/t — with 
n = 0.75, 0.83, 0.88 and 0.92. (Right panel) t'/t = -0.1 with 
n = 0.80, 0.83, 0.85 and 0.87. The position of the van Hove 
singularity (peak in the density of states) shifts from positive 
frequency to negative frequency as the filling moves towards 
half-filling. The quantum critical point is at n c ~ 0.88 for 
t'/t = 0, and n c w 0.835 for t'/t = -0.1. 



also been observed in the "momentum-selective" metal- 
insulator transition scenario^—, however, quantum crit- 
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ical phenomena associated with the van Hove singularity 
and the Lifshitz transition have not been discussed there. 
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FIG. 3. (color online) Intensity plots of the spectral function 
A(k,ou) along the segment M — ► X —> F of the symmetric 
path inside the Brillouin zone. For different t' /t near the 
quantum critical fillings: t' ft = with n = 0.88 (left), t' ft = 
-0.1 with n = 0.83 (middle), and t'/t = -0.2 with n = 
0.78 (right), a saddle point region (flat dispersion) around X 
crosses the Fermi level at the critical filling n c . The region 
becomes wider as t' goes to negative values. 

The van Hove singularity originates from the flat dis- 
persion along the antinodal direction with energy close to 
the Fermi leve l 36 i 47 i 48 , as shown in Figure [3l At the criti- 
cal dopings for t' jt = 0.0, —0.1, and —0.2, there is always 
a saddle point region around the momentum X = (7r, 0) 
where the dispersion is flat. In addition, we find that 
the flat region become wider and more pronounced as 
t' becomes negative. Such flat dispersion at the chem- 
ical potential contributes low energy states and results 
in the van Hove singularity at the Fermi energy. Hence, 
our observation of the Lifshitz transition at the quantum 
critical doping at t'/t < is closely tied (if not in one- 
to-one correspondence) to the crossing of the Fermi level 
by the van Hove singularity. 

The next-nearest-neighbor hopping t' can be viewed as 
the control parameter of the Lifshitz transition. To em- 
phasize this point, we also calculate the zero frequency 
spectra for the non-interacting system as a function of 
t' It, and compare them with our DC A results for the 
interacting system. In Figure 2J the left panel demon- 
strates the effect of t' on the Fermi surface topology of 
the non-interacting system. We find the Fermi surface 
changes from electron-like to hole-like when t' goes from 
positive to negative values. The right panel of Fig. [H 
collects our results for the spectra of the interacting sys- 
tem, with n = 0.85 fixed and t' /t varies from 0.3 to -0.2. 
Similar to the noninteracting system, the Fermi surface 
topology also changes from electron-like to hole-like as 
t' jt goes from positive to negative, which confirms that 
t' is indeed the relevant control parameter of the Lifshitz 



transition in the interacting system. However, the inter- 
acting system follows a more complicated phenomenol- 
ogy than the non-interacting system, due to the effect of 
electron-electron interaction. For example, in the lower 
two panels of the right hand side of Fig. 2J where t' is 
negative, not only is the Fermi surface topology changed 
to hole-like, but the spectral weight along the antinodal 
direction (7r, 0) also vanishes, signaling that the system 
enters the pseudogap phase. It is the strong interaction 
which results in a redistribution of the spectral weight 
and leading to features different from the non-interacting 
system. 

Instead of obtaining the spectral information from the 
analytical continued data, one can also directly read off 
the quasiparticle weight Z(k) from the Matsubara fre- 
quency data. Since the quasiparticle weight will be fi- 
nite across a Fermi surface, but vanishes if the spectrum 
is incoherent, it can be used to distinguish a Fermi liq- 
uid and a pseudogap state. The quasiparticle weight is 
calculated from the Matsubara frequency self-energy as 

S"(fc,ia;o) \ — 1 1 

I k—k f ' 



Z (k F ) = (1 



where ujq = irT is 



the lowest fermionic Matsubara frequency. In the limit 
T — >> 0, and for the Fermi- liquid self-energy, Zo(fe) con- 
verges to the quasiparticle weight, Z(k). Figure [5] shows 
Zan(T) = Zo(uo = 7rT, kp), the Matsubara quasiparti- 
cle weight on the Fermi surface along the antinodal direc- 
tion for various t' jt. Zan exhibits different behavior for 
n < n c where the quasiparticle weight approaches a finite 
value and for n > n c where the quasiparticle weight van- 
ishes in the limit T — >• 0. The temperature dependence 
of Zan(T) furthermore provides information about the 
relevant energy scales. One sees a clear difference in the 
temperature dependence of Zan(T) when the system en- 
ters Fermi liquid, marginal Fermi liquid and pseudogap 
phases. Inside the Fermi liquid phase, n < n c , where we 
choose n = 0.75 for t' jt = 0, n = 0.70 for t' jt = -0.1 
and n = 0.65 for t'/t = —0.2, Zaf(T) extrapolates to a 
finite value roughly linearly at low T. Close to the Lif- 
shitz transition points, n « n c , where we choose n = 0.88 
for t'/t = 0, n = 0.83 for t'/t = -0.1 and n = 0.78 for 
t'/t = —0.2, Zan(T) shows a behavior consistent with 
the marginal Fermi liquid picture, i.e., a negative curva- 
ture in Zaf(T) at low T— . Finally, when the system is 
inside the pseudogap phase, n > n c , we choose n = 0.95 
for t' = 0, n = 0.95 for t'/t = -0.1 and n = 0.95 for 
t'/t = —0.2, Zaf(T) goes to zero even faster than that 
in the marginal Fermi liquid region. The detailed in- 
formation about the crossover temperatures T* and Tx , 
where the former signifies the crossover between pseudo- 
gap and the marginal Fermi liquid and the latter is the 
crossover temperature between the marginal Fermi liquid 
and the Fermi liquid, has been discussed in detail in our 
previous publication s 20 ! 23 . 

The bare d-wave pairing susceptibility is calculated 
as Xod(u) = EfcXo(^,4 = 0)g d (k) 2 /^2 k g d (k) 2 , where 
9d(k) = (cos(k x ) — cos(ky)) is the d-wave form factor. 
It exhibits significantly different features near the QCPs 
compared with those predicted by static van Hove sin- 
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FIG. 4. (color online) t' jt as the control parameter of the Lifshitz transition. (Left panel) A(k,oj — 0) for the non-interacting 
system at half- filling. As t' /t goes from positive to negative, the Fermi surface changes from electron- like to hole- like. (Right 
panel) A(k,uj = 0) for the interacting system, U/t — 6, with filling n — 0.85 fixed. At positive t\ the system is inside the 
metallic Mott gas phase, and has complete Fermi surface; at negative t', the system is inside the pseudogap phase, and the Fermi 
surface becomes hole- like. Since the Coulomb interaction reshuffles the spectral weight, at t' /t = —0.1, —0.2, the pseudogap at 
antinodal direction can be clearly seen. 



gularitie o 35 ! 36 . As shown in Fig. [6j close to the quan- 
tum critical fillings, the real part of the bare pairing sus- 
ceptibility x'od^ 00 — 0,T) diverges quickly with decreas- 
ing temperature, following a power-law behavior close 
to 1/y/T, which is different from the BCS type of log- 
arithmic divergence Xq(T) oc N (0) \ii(uj jj / T) with N(0) 
the single-particle density of states at the Fermi surface 
and ujd the phonon cutoff frequency. Such behavior is 
consistent with our previous results at t' = close to 
the quantum critical fillin g 22 ! 24 . More interestingly, we 
find at t' < 0, such temperature dependence persists. 
As shown in Fig. [6j the prefactor, a, in the power- 
law term, becomes larger as t' becomes more negative 
. It signifies that the divergence becomes stronger at 
t' < 0, and if the pairing strength does not change, there 
will be a higher superconducting transition temperature 
T c compared with that in the t' = case. Moreover, 
when we scale the imaginary part of the bare pairing 
susceptibility as T 1,5 x"od(kO/^ versus cj/T 49 , as shown 
in Fig. [7J we find the curves from different tempera- 
tures fall on the same universal scaling function such that 
T L VodM/^ = H{u)/T) * (u/T)- 1 - 5 for uj/T > 9 « 
4t/J, where J w Q.l l 20 i 50 is the antiferromagnetic ex- 
change energy near half filling. From the Kramers-Kronig 



relation, the real and imaginary parts of the susceptibility 
are related via Xo^(^) = ^ f < ^ Co 'X"o(i(^)/^, so collapse of 
the x"od(^)/^ will contributes a term oc -j= in the real 
part of the Xod( T )- 



IV. DISCUSSION 

One of the important issues in the field of cuprate 
physics is the identification of the quantum critical point 
accepted by and large as relevant for the superconduc- 
tivity and in particular the strange "non-Fermi liquid" 
behavior observed in the vicinity of the maximum of 
the superconducting dome. Based on our numerical re- 
sults presented here together with those published previ- 
ously 2 ^ - —, we suggest the schematic phase diagram for 
the two dimensional Hubbard model shown in Figure [8] 
(note that for clarity we purposely neglect the super- 
conducting and the antiferromagnetic transitions in this 
schematic plot). The control parameters are the next- 
nearest-neighbor hopping t' It, chemical potential /i, and 
temperature T. For each t' ', we shift the chemical \i with 
respect to its value at half filling fiuF(t f /t). For t' > 0, 
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FIG. 5. (color online) Matsubara quasiparticle weight 
Zan(T) versus temperature T evaluated with k on the Fermi 
surface along the antinodal direction, for various fillings and 
t' . For n < n c , Zan(T) displays the Fermi liquid behavior, 
with Zan(T) extrapolating linearly to a finite value at low 
T; for n ~ n c , Zan(T) displays the marginal Fermi liquid 
behavior, with negative curvature in Zan(T) at low T; and 
once inside the pseudogap region, n > n c , Zan(T) goes to 
zero faster than that in the marginal Fermi liquid fillings as 
a function of T. 



FIG. 6. (color online) The real part of the bare <i-wave pairing 
susceptibility, Xodi 00 = 0,T), at zero frequency, close to the 
quantum critical fillings for t'/t = 0, t'/t = —0.1, and t f /t = 
-0.2. The solid lines are fits to x'od(u = 0,T) = a/VT + 
b log (c/T). Close to the quantum critical fillings (n = 0.85 
for t'/t = 0, n = 0.83 for t'/t = -0.1, and n = 0.78 for 
t'/t = —0.2), x'odi 00 — 0, T) shows power-law divergence with 
decreasing temperature. The prefactor a associated with the 
square-root term increases as t' become negative, signifies a 
stronger divergence in the pairing susceptibility. 



there is a first order phase separation transition between 
an insulating and incompressible Mott liquid (ML region 
in the lower inset of Fig. [8j) and a weakly compressible 
metallic Mott gas (MG region in the lower inset of Fig. [8]). 
As a function of temperature, this region of first order 
transitions terminates in a line of second order classi- 
cal critical points (dots and line for t' > 0). By vary- 
ing doping and t', this line of second order transitions is 
suppressed to T = 0, the corresponding T = 0, t' = 
doping, defines the quantum critical point in the phase 
diagram. This QCP separates the pseudogap regime from 
the Fermi liquid. As usual, one observes a V-shaped 
quantum critical region above this QCP, which shows the 
signatures of the marginal Fermi liquid. Note that this 
QCP at t f = also involves a Lifshitz transition, i.e. a 
change in Fermi surface topology from electron-like in the 
Fermi liquid phase to hole-like in the pseudogap phase. 
This property can be directly inferred from the changes 
in the low-frequency structures of the spectral function 
(see Fig.[T]). At the critical doping the van Hove singular- 
ity in the single particle density of states is located at the 
Fermi energy (see Fig. [2]). Note that there are actually 
two effects involved: Firstly, the change in the overall 
shape of the Fermi surface, and secondly a distinctive 
suppression of spectral weight along the antinodal direc- 
tions in the pseudogap phase. Based on the momentum 
resolution presently available, we are not able to make 
decisive statements about a possible formation of hole 
pockets. This last step would be necessary to unam- 
biguously establish the nature of the Lifshitz transition. 
We are currently developing a multi-scale dual Fermion 
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FIG. 7. (color online) Plots of T^V'odM/cJ versus u/T 
close to the quantum critical dopings (n = 0.85 for t'/t = 0, 
n = 0.83 for t'/t = -0.1, and n = 0.78 for t'/t = -0.2). As 
the temperature decreases, the curves coincide for uo/T > 9 ~ 
(At/ J) defining a scaling function H(lj/T), corresponding to 
a contribution to Xod(T) = \ f dux" od(u) / uo oc A= as shown 
in Fig. El In the scaling regime, H(uj/T) & (cj/T)" 15 . 

dynamical cluster approach^, which combines the DCA 
method used in this work with the recently introduced 
dual-fermion formalism 52 . This approach can systemat- 
ically incorporate the long-ranged correlations through 
the dual fermion lattice calculation, and will eventually 
provide us the necessary momentum resolution to address 
the question of the hole pockets and the type of the Lif- 
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FIG. 8. (color online) Schematic phase diagram of the 
Hubbard model close to the quantum critical point (QCP) 
with temperature (T), chemical potential (/i), and the next- 
nearest-neighbor hopping (£') as the control parameters. For 
clarity we neglect superconducting and antiferromagnetic 
phases and only focus on the Lifshitz physics. For each t', we 
shift the chemical potential fi with respect to its half filling 
value /iHF^/t). The three insets show the cut of the phase 
diagram at three different t' /t = 0.1, —0.1 and —0.2. At posi- 
tive t\ the Mott liquid and Mott gas phases are separated by 
a first order line; at negative t\ the pseudogap and the Fermi 
liquid phases are separated by a quantum critical region, and 
the region becomes wider as t' moves towards the negative di- 
rection. As the Lifshitz points move to higher hole doping at 
negative t\ the slope of the boundary between the pseudogap 
and the quantum critical region becomes less steep. 



shitz transition associated with the QCP. However, the 
numerical data presented in this work are very indicative 
and indeed consistent with a collapse of the Fermi surface 
from a large to a small one. 

As t' becomes negative, phase separation is completely 
suppressed, but the critical behavior associated with the 
Lifshitz transition remains, leading to a line of critical 
points at T = 0, as sketched in the schematic phase dia- 
gram in the region t' < 0. Since the effect of negative t' 
is to distort the Fermi surface towards hole-like (see Fig. 
2]), in the t 1 jt < part of the phase diagram, one nec- 
essarily needs higher hole doping concentration to have 
the van Hove singularity locate at the Fermi level. This 
provides us a simple explanation for the fact that the 
Lifshitz points move to higher hole doping at t' < 0. 
Across the Lifshitz points, it is not only the change in 
shape of the Fermi surface, but more importantly the 
strong suppression of weight along the antinodal direc- 
tion which leads to the interpretation of a Lifshitz tran- 
sition at T = 0. The quantum fluctuations here again 
lead to a region with quantum critical behavior. As is 
commonly accepted, the low-energy model for the high- 
T c cuprates is the two dimensional Hubbard model with 
at least t > and t' < 0, we thus have strong numerical 



evidence that the QCP is indeed due to a transition of the 
Fermi surface topology, as already suggested by several 
other authors 1 "—. We also find at negative t' the quan- 
tum critical region becomes wider in the doping range 
(see the two upper insets of Fig. [8]), the quantum critical 
points move to higher hole doping concentrations, and 
the boundary between the pseudogap region and quan- 
tum critical region becomes less steep as t' moves towards 
the negative direction. 

There is, however, an important missing link: The ex- 
planation, for the interacting electron system, of why the 
van Hove singularity crosses the Fermi surface at the crit- 
ical filling n c , even for t' = 0, thus triggering the Lifshitz 
transition. In previous DCA studies of the 2D Hubbard 
mode l 3Q i 53 , it was found that the non-local correlations 
present in the DCA, but missing the DMFA, tend to dis- 
tort the Fermi surface of the hole doped system so that 
it is centered around the wavevector (7r,7r) rather than 
(0,0). As the system is doped away from half filling, 
these correlations become weaker, and the DCA Fermi 
surface returns to the DMFA (or uncorrelated) Fermi 
surface centered around (0,0). Thus the topology of the 
Fermi surface changes from hole-like to electron-like at 
finite doping. 

On the other hand, since the correlation effects at finite 
doping show strong temperature dependence, we expect 
the phenomena associated with the van Hove singular- 
ity also have a temperature dependence which is quite 
different from those predicted by models with a static 
van Hove singularity resulting from saddle point in the 
bare dispersion 3 ^. This has important consequences for 
the superconducting transition T c , which, at sufficiently 
high doping, is determined by the BCS-like condition 

VdXodi 00 = 0) = 1 where Xod * s ^ ne rea ^ P ar ^ °f ^ ne bare 
d-wave pairing susceptibility and Vd is the strength of the 
d-wave pairing interaction. 

Following the results in Figure [6] and Figure [7J we find 
at t' = 0, t'/t = -0.1, and t'/t = -0.2, x'odi^ = °) 
close to the van Hove singularity diverges quickly with 
decreasing temperature, roughly following the power-law 
behavior ex a /y/T, which is significantly different from 
the conventional BCS logarithmic divergence or the log- 
squared divergence found at a static van Hove singular- 
it y 35 i 36 . Moreover, we find that the prefactor a increases 
with negative t' \ which signifies a stronger divergence in 
the pairing bubble compared with the t' = case. How- 
ever, we also know the d-wave pairing vertex, Vd, de- 
cays monotonically with hole doping, since Vd originates 
predominantly from the spin channel 22,54 . A negative 
t' will frustrate the antiferromagnetic background, and 
hence suppress the spin channel contribution to the d- 
wave pairing vertex. At the same time, according to 
our Lifshitz phase diagram (Fig. [8]), one needs higher hole 
doping concentration to approach the QCP for negative 
t' . The higher quantum critical doping and the suppres- 
sion of the spin channel contribution to the pairing ver- 
tex at critical doping lead us to expect the d-wave pairing 
vertex, Vd, becomes even smaller for t f < than it's t f = 
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correspondence. Hence, a more strongly divergent pair- 
ing bubble, x'odi anc ^ a wea ker pairing vertex, Vd, bring 
us to an interesting situation for t' < 0, where there is a 
competition between x'od anc ^ Whether this compe- 
tition will yield a higher T c and larger superconducting 
dome in the doping range close to the Lifshitz points at 
negative t' will be the subject of our future investigation. 

Our findings of the phase diagram are also phenomeno- 
logically consistent with the recent proposal of unconven- 
tional quantum criticality which appears at the border of 
the first order metal-insulator transition and continuous 
transition triggered by quantum fluctuations^—. This 
border, the marginal quantum critical point, belongs to 
an unprecedented universality class and has unique fea- 
ture with combined characters of symmetry breaking and 
topological(Lifshitz) transitions^. Whether the Lifshitz 
points presented at t'/t < in this work obtain uncon- 
ventional quantum criticality is also a subject for our 
future study. 

V. CONCLUSION 

Using large-scale dynamical cluster quantum Monte 
Carlo simulations, we map out the Lifshitz phase diagram 
of the two dimensional Hubbard model in the vicinity of 
the quantum critical fillings. The control parameters of 
the phase diagram are temperature T, chemical potential, 
and the next-nearest-neighbor hopping t f /t. Consistent 
with our previous result a 21 i 23 i 24 , we find at the positive 
t', there is a first order phase separation transition which 
is terminated by a second order critical point. As t' — >• 0, 
the second order terminus is driven to zero temperature, 
and becomes the QCP separating the pseudogap phase 
and the Fermi liquid phase. Here, we extend the inves- 
tigation into negative t' , and find out for t'/t < 0, at 
the ground state, there is a line of van Hove singulari- 
ties, at which the Fermi surface topology changes from 
electron-like in the Fermi liquid region to hole-like in the 
pseudogap region of the phase diagram. The points on 
this line of van Hove singularities hence are the quan- 



tum critical points where the Lifshitz transition occurs. 
Close to these QCPs, the bare d-wave pairing bubble 
diverges algebraically in temperature. Originating from 
these Lifshitz points, the V-shaped quantum critical re- 
gion emerges with the marginal Fermi liquid properties 
and vanishing quasiparticle weight 20,24 . We also find the 
V-shaped quantum critical region becomes wider in dop- 
ing range (and chemical potential) as t'/t becomes nega- 
tive. 

There remain a number of interesting open issues, in- 
cluding the possible formation of hole pockets in the pseu- 
dogap phase, the explanation of the van Hove singularity 
crossing the Fermi surface at the quantum critical dop- 
ing, the superconducting transition temperature T c and 
the shape of the superconducting dome at negative t', 
and the possible unconventional quantum criticality as- 
sociated with the Lifshitz transition points. All these 
questions require not only massively parallel large-scale 
simulations, but more importantly, new techniques that 
will greatly increase our momentum resolution. In fact, 
progress has already been made along these directions^ 1 -, 
and the remaining questions will be addressed in future 
work. 
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